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Abstract. We study the effects of vertical outflows on mean-field dynamos 
in disks. These outflows could be due to thermal winds or magnetic buoy- 
ancy. We analyse numerical solutions of the nonlinear mean-field dynamo 
equations using a two-dimensional finite-difference model. Contrary to ex- 
pectations, a modest vertical velocity can enhance dynamo action. This 
can lead to super-exponential growth of the magnetic field and to higher 
magnetic energies at saturation in the nonlinear regime. 



1. Introduction 

Large scale magnetic fields are often considered to be an important factor 
in accretion disks (including generation of turbulence) , crucial in launching 
winds or jets. The origin of such large scale magnetic fields in accretion 
disks is still unclear: they may be advected from the surrounding medium 
or be generated by a dynamo inside the disk. Advection appears unlikely 
as turbulence leads to enhanced viscosity and magnetic diffusivity, so that 
the two are of the same order, i.e. the magnetic Prandtl number is of order 
unity [1]. In this case the turbulent magnetic diffusivity can compensate the 
dragging of the field by viscously induced accretion flow [2, 3, 4]. Turbulent 
dynamo action is a plausible mechanism for producing large scale magnetic 
fields in accretion disks [5, 6]. Dynamo magnetic fields can launch winds 
from accretion disks [7, 8, 9]. However, the wind can also affect the dynamo. 
In particular, a dynamo enhancement by winds was suggested earlier for 
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galactic dynamos [10]. In the same context the effects of shear in the vertical 
velocity on the dynamo was considered [11]. 

We assume here that the magnetic fields are generated by a dynamo 
acting in a relatively thin accretion disk. We show how vertical flow can 
enhance the dynamo, allowing for a larger growth rate and leading to super- 
exponential growth of the magnetic field and to enhanced saturation levels 
of magnetic energy. 

The vertical velocities can have several origins. For example, they can 
be due to a thermally driven wind emanating from the disk or magnetic 
buoyancy in the disk. Note that here we invoke magnetic buoyancy as a 
driver of vertical outflows; it can itself contribute to the dynamo effect [12], 
but such effects are not taken into account here. 

2. The model 

The equation which we solve is the mean-field induction equation which we 
evolve in terms of the vector potential A, where B = V x A, 



where j = V x B//xo is the current density, no the magnetic permeability, r\ 
the turbulent magnetic diffusivity, a the a-effect and V the mean velocity. 
We neglect the radial component of the mean velocity. Assuming the a- 
and 77-tensors to be isotropic, we can consider scalar quantities. We do 
not make the thin-disk approximation but solve the general equations. We 
adopt cylindrical coordinates (r, <p, z) and restrict ourselves to axisymmetric 
solutions. Equation (1) is solved using a sixth order finite-difference scheme 
in space and a third order Runge-Kutta time advance scheme. 

Our computational domain contains a disk embedded in its surrounding 
halo. We take a disk aspect ratio of /idisk/-Rdisk = 0.1 and a halo with 

^halo/^disk ~ 6. 

As an initial condition for the magnetic field we choose a purely poloidal 
field in the disk of either dipolar or quadrupolar symmetry. 

On the boundaries of the computational domain we impose pseudo- 
vacuum conditions. However, since the boundaries are in the halo far away 
from the disk, the choice of boundary conditions is not crucial. 

The a-coefficient a(r, z) is antisymmetric about the disk midplane and 
vanishes outside the disk. We adopt 



— - = Vx(VxA)+aVxA- Wo j, 



(1) 




for \z\ < h, 
for \z\ > h. 



(2) 



The ^-profile cuts off the a-effect at the outer radius of the disk as well as 
at an inner radius, where the rotational shear is very strong. 
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As appropriate for accretion disks, we adopt a softened Keplerian an- 
gular velocity profile in r in the disk as well as the halo, 

where G is Newton's gravitational constant, M is the mass of the central 
object, ro = 0.05 is the softening radius, and n = 5. At r = 0, VL vanishes. 
The turbulent magnetic diffusivity is given by 

7](r, z) = 7] hsdo + (r/ disk - r/ halo )£(r, z). (4) 

The profile £(r, z) defines the disk: £ is equal to unity inside the disk and 
vanishes in the halo but with a smooth transition between. We carried out 
computations for two cases: homogeneous conductivity, Tyhaio/^disk = 1> an d 
low conductivity in the halo, r/haio/^disk = ^0- 

The profile for V z is linear in z, V z (z) = V z oz/h, where V z $ is a charac- 
teristic vertical velocity. 

Our dynamo problem is controlled by the three magnetic Reynolds num- 
bers related to the a-effect, the differential rotation and the vertical velocity, 

R a = ha /r] disk , R u = h 2 S/rj disk , R v = hV z0 /r} Aisk , (5) 

where S(r) = rdQ/dr « —3/2 Q(r) is the rotational shear. The dynamo 
number is defined as T> = RaR^. Since the dynamo number is approxi- 
mately constant, \V\ ~ agg ~ fO 2 — f0 4 in thin accretion disks [5], we 
assume V to be constant by setting a>o(r) = Vr]^ isk / S(r)h 3 . 

Because of the strong differential rotation we assume that the magnetic 
field in the disk is generated by a standard aoj-dynamo. The two con- 
trol parameters are then the dynamo number T> and the vertical magnetic 
Reynolds number R v . The value of R v obviously depends on the origin 
of the vertical outflow. With magnetic buoyancy, a rough estimate can be 
done assuming that the vertical velocities are comparable to the Alfven 
speed. Estimating the Alfven speed from magnetic equipartition, one gets 
a vertical magnetic Reynolds number R v of order unity. 

3. Linear results for R v = 

F° r ^halo = f/disk) the growth rates of dipolar and quadrupolar modes should 
be exchanged when T> reverses sign, i.e. the graph of Re(-y) as a function 
of V should be symmetric with respect to the vertical axis V = with 
quadrupolar and dipolar modes exchanging their roles [13]. Our numerical 
simulation shows this symmetry with a good precision (Fig. la). This is 
because the ratio between /ihaio and fadisk is large enough, about 6. 
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Figure 1. Real part of the growth rate of the magnetic field as a function of —T> 
with R v — 0. Asterisks denote non-oscillatory, diamonds oscillatory solutions. Solid lines 
denote quadrupolar (Q), dashed lines dipolar (D) modes, (a) is with homogeneous con- 
ductivity, (b) for low-conducting halo. Note the symmetry of (a) with respect to V = 0. 



For V < (R a > 0) the first leading growing mode is quadrupolar 
whereas for V > (R a < 0) it is dipolar; both are non-oscillatory. We 
thus study the effect of vertical velocities for T> < in quadrupolar and for 
V > in dipolar symmetry at dynamo numbers up to ±300. 

For rj = const., the critical dynamo number is | Dent I — 5. The dominant 
mode changes its symmetry and becomes oscillatory at \T>\ ~ 300. 

The qualitative behaviour in the case of low conductivity in the halo 
is quite similar to that of homogeneous conductivity. Also the quadrupolar 
non-oscillatory modes first become dominant for R a > and the dipolar 
non-oscillatory modes for R a < (compare Figs, la and b). 

But there are quantitative differences. The value of \T>\, where the sym- 
metry of the leading mode changes, is larger than 700. Also, the diagram 
in Fig. lb is not symmetric. The critical dynamo number is V ~ —50 (50) 
for quadrupolar (dipolar) modes. 

4. Linear behaviour with vertical velocities 

As discussed in section 3, we only consider quadrupolar symmetry for V < 
0. As shown in Fig. 2, the growth rate of the magnetic field Re(7) is a 
non-monotonous function of R v , for ??halo/^disk equal to both 1 and 20. A 
maximum in Re(7) occurs, however, only for \T>\ large enough, and the 
smaller the magnetic diffusivity of the halo, the smaller is the required 
\T>\. The larger |£>|, the more pronounced is the maximum. The maximum 
growth rate occurs for R v of the order of 1 to 10. 

In Fig. 3 we show the effect of the dynamo number and the vertical 
velocity on the magnetic field configuration for r/halo = ?7disk- Increasing \V\ 




from 50 to 200 in the absence of any vertical velocity results in the magnetic 
field becoming concentrated at larger radii in the disk (Fig. 3a and b). 

As can be seen in Fig. 3b, c and d, a vertical velocity with R v < 7.5 leads 
to a wider vertical distribution of magnetic field which reduces magnetic 
diffusion and enhances the dynamo action. At still larger R v , the wind 
aligns the poloidal field lines with the lines of constant rotation, and the 
fi-effect is slowly switched off leading to a decrease in the growth rate, 
which eventually becomes negative. The mode structures at a maximum of 
Re(7), reached at R v = 7.5 for V = —200) and R v = 3 for V = —50, are 
very similar to each other. 

The growth rate for dipolar non-oscillatory modes (V > 0) has no max- 
imum as R v increases; at R v of order 1 or even less the dynamo is switched 
off. The magnetic field is advected outwards and aligned with the O-contour 
lines in vertical direction very quickly. 

5. Nonlinear behaviour with vertical velocities 

We consider nonlinear, saturated solutions in a model with homogeneous 
conductivity, r/h a io = f?disk, negative dynamo number and quadrupolar sym- 
metry (the dominant symmetry for moderate \T>\ with R v = 0). All satu- 
rated magnetic fields are non-oscillatory, as the corresponding linear modes. 

5.1. MAGNETIC BUOYANCY 

We parameterize the effect of magnetic buoyancy by assuming the vertical 
velocity to be proportional to the maximum magnetic field strength. Thus, 
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Figure 3. The effect of dynamo number and vertical velocity on the magnetic field 
configuration, for the linear model with ?7haio = ?7disk. Shown are the poloidal field lines 
(solid is clockwise, dotted is counter-clockwise) and the toroidal field (grey scales: bright is 
in positive, dark in negative azimuthal direction), (a): T> = —50, R v = 0; (b): T> = —200, 
R v = 0; (c): V = -200, R v = 0.75; (d): V = -200, R v = 7.5 (where Re(7) is maximum). 
Comparing (a) and (b) shows the effect of increasing T> while R v = 0. Panels (b) to (d) 
follow the upper curve in Fig. 2a. All modes are quadrupolar and non-oscillatory. 



the vertical magnetic Reynolds number becomes time-dependent, 

R v = R v omax(B r ,B v )/\B \, R v0 = V z0 h/r] disk , (6) 

where Bq is a characteristic field strength. 

As time increases, |B| grows and therefore R v is increasing and the 
growth rate follows, e.g., the upper curve in Fig. 2a for T> = —200. Thus, 
as long as R v is less than the position of the 7-maximum, 7 increases with 
time, which results in a super-exponential growth. After the maximum, 
at R v = 7.5, 7 decreases and eventually the magnetic field approaches its 
saturation level. This is reached when R v = R m , where R m is the zero of 
7(i£u). According to Eq. (6), the saturation value of magnetic energy will 
thus be £ mag cx 1/R% . 
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Figure 4- Magnetic energy as a function of R v for the model with a-quenching (left) 
and for the models with a-quenching and/or magnetic buoyancy (right) for ??haio = f7disk- 
All modes are quadrupolar and non-oscillatory. The right panel shows the case V = —50. 



5.2. a-QUENCHING 

We consider the back-reaction of the magnetic field on the a-effect, intro- 
ducing the nonlinearity 1/(1 + B 2 /B^ q ) as a factor in the a-effect; B e<i is 
the equipartition field. 

The magnetic energy has a maximum at a certain value of R v (Fig. 4, 
left), but this maximum occurs at R v ~ 0.1, a value smaller than where 
a maximum of Re(y) occurs in the linear regime, and it is also less pro- 
nounced. The value of R v where the maximum occurs, is independent of V. 
The magnetic energy scales roughly with the dynamo number, E mSug oc T>. 

The effect of vertical velocity on the magnetic field for the model with 
a-quenching is shown in Fig. 5; the structure changes only weakly with 
the dynamo number. With increasing R v the poloidal field lines become 
vertical, i.e. aligned with the lines of constant rotation. Hence, shear has 
no effect, and the fi-effect is switched off, as in the linear model. 

5.3. a-QUENCHING TOGETHER WITH MAGNETIC BUOYANCY 

When a-quenching and magnetic buoyancy are combined, the vertical mag- 
netic Reynolds number depends on position and time and takes the form 
R v = R v q^J B 2 + B 2 /\B cq \. Without a-quenching, since R v is now nonlocal, 
the temporal behaviour of the growth rates at short times does not have to 
be super-exponential. Surprisingly, again we find E mag oc 1/R% (see Fig. 4, 
right, dashed line). 

In the second panel of Fig. 4 magnetic energy is plotted against R v for 
V = —50 for the three nonlinear models, a-quenching appears to be the 
dominant nonlinearity in the model considered. 
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Figure 5. The effect of vertical velocity on the magnetic field configuration for the model 
with a-quenching and ryhaio = ?7disk- Shown are the poloidal field lines and the toroidal 
field, (a): V = -50, R v = 0; (b): V = -50, R v = 3; (c): V = -50, R v = 7.5. Panels (a) to 
(c) follow the lower curve in Fig. 4, left. All modes are quadrupolar and non-oscillatory. 
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